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1 Introduction 



Many problems of magnetohydrodynamics do not have analytic solutions and have to be ap- 
proached using simple models that help to understand some important phenomena and to distinguish 
amongst the more and less influential ones. 

Terrestrial magnetosphere is the layer of the atmosphere related to terrestrial magnetic field, 
where interactions between the magnetic field and the wind coming from the Sun occur, deforming 
the magnetic field in the direction of solar wind flow. The solar wind is a " ionized gas" with a very 
large electrical conductivity, chiefly made up of protons and electrons of high velocities ^ . 

As the velocity of this fluid is so large, far away from the source of magnetic field it is not affected 
appreciable due to terrestrial magnetic field, however the magnetic field changes the trajectory of 
particles near the Earth, partly this is the cause that certain low energy particles remain trapped 
in the magnetic field oscillating and producing electromagnetic radiation, and in the poles their 
interaction with air molecules cause the know effect called northern dawn 

This paper is organized as follows: In section 2 we derivate the Magnetohydrodynamic equations, 
in section 3 we describe the problem and calculate the deformation of terrestrial magnetic field due 
to the solar wind. Section 4 is devoted to analysis of the results and conclusions. 

2 Magnetohydrodynamic 

The laws of magnetohydrodynamics j^j are deduced from Maxwell's equations: 

V-i? = 47r— . (2. f) 



IdB , , 

Vxi..^-- (2.2) 

V • -B = (2. 3) 

WxB^—fij + ^— (2.4) 

Introducing the Ohm's law for a conductive fluid, i.e. that the electric current circulating inside the 
fluid is proportional to the electric field, this is: 

f - aE' (2. 5) 

where j' and E' are the electric current and the electric field as seen from the frame of reference of 
the fiuid. By mean the Lorentz's transformations we pass from fluid reference frame to laboratory 
reference frame, in the non-relativistic approximation w << c. 

J = /' + PeV (2. 6) 



E' ^ E +- y. B (2. 7) 

c 
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Assuming a conductive fluid without net charge, i.e., the charge density pe is nufl. By substituting 
equations p. 6|l and p. 7|l into equation p. 5() we obtain the final expression for the electric current: 



j = a[E 



X B 



If we replace equation (|2. 8|l into equation (|2. 4() . we obtain: 



47r 

y X B — — 

c 



E+ - X B 
c 



ep, dE 



(2. 



(2. 9) 



In this equation we observe that the second term on the righthand side can be despised, since its order 
of magnitude is very small in comparison to the first term as the fluid is a conductor (conductivity 
a very large). With this consideration equation (|2. 9|l becomes: 



47r 

V X B ~ — /i(T 

c 



E + - X B 
c 



(2. 10) 



We apply the rotational operator to equation 1)2. lUI) and, keeping in mind that the divergence of 
the field is always null, we obtain: 



, - An 
V^B = —fia 
c 



X B 



(2. 11) 



V X £; + V X 

Consider Faraday's law equation H2. 2|) we arrive to a diffusion equation for the magnetic field: 

(2. 12) 



dB , - , 

— = Tj^y^B + \/x{vxB) 



where rjri 



is known as the magnetic viscosity. The diffusion equation must be accompanied 



with the movement equation (Newton's second law) for the fluid and the equation of continuity [S]: 

(2. 13) 



dv 1 ^ ^ o 

p— = -VP + -j X B + yyV^w. 
dt c 



V • ipv) = (2. 14) 

Equation 12. 141) is the continuity equation for mass, where p is the density of matter and v the 
velocity of the fluid; equation H2. 13() is called the Navier-Stokes's equation or movement equation for 
fluids. In this equation intervenes the pressure gradient, a term of magnetic force and a term due to 
the viscosity of the fluid. If we considered that the fluid is incompressible, equations (|2. 12|I - H2. 14|) 
are a complete set describing the dynamics of the fluid and the magnetic field, on the other hand, 
if the fiuid is considered compressible, in addition of equations H2. 12|) - (|2. 14|l the state equation of 
the fiuid should be in mind. Equations H2. 12() - (|2. 14|) are the magnetohydrodynamic equations. 



3 Problem Description 

The objective of this work is to analyze how the terrestrial magnetic field is modified by the 
presence of the solar wind in zones far away from Earth's magnetic field source. For this, we consider 



3 



a conductive fluid whose velocity is constant, i.e. v — vqCLx. This approximation is valid in the sense 
that the velocity of solar wind particles is very high far away of Earth, and is not affected appreciably 
by the terrestrial magnetic field. Besides, we also consider an incompresible and irrotacional fluid, 
which is equivalent to exclude the dynamics related to the formation of vortices in the fluid, that 
generally form in the proximities of Earth, because its magnetic field is very strong. According to 
these considerations, the velocity of the fluid can be written as the gradient of a potential function 
of the form: 

v = -V(j)y. (3. 15) 

From equation p. 14() . for an incompressible fluid its density is constant, then we obtain: 

V • (pv) = 

p(V • w) + Vp • t/ = 

V • V = 

V^0„ = 0. (3. 16) 

Another aspect that we consider is the stationary case in which the physical quantities remain 
constant in time. This approximation can be achieved if we take a constant velocity for the solar 
wind. Another factor might carry changes in Earth's magnetosphere is related to possible variations 
of Earth's magnetic dipole moment, however it has a very small effect during small intervals of time 
0]. Accordingly, equation (|2. 2|l for the electric field can be expressed as: 

V X E = 0. (3. 17) 

In this model we also suppose that the fluid does not possess net electrical charge as, although the 
gas in consideration is composed by charged particles, to great scale the net charge is very small. 
Then Gauss's Law ec. (|2. l|l is given by: 

V -E^O. (3. 18) 
Taking into equation (|2. ll|l V x £^ = 0, we obtain 

„9 ^ 47r I vx b\ 

this last expression is known as the diffusion equation for magnetic fields, which can be obtained 
also by considering that the electric field is null, i.e., if we take a the trivial solution for equation 
(j3. 18|l . then we have E = 0. In our model obtaining a non-trivial solution for the diffusion equation 
(|2. 12|l can be a complicated calculation, as it also has to verify the condition W ■ B — Q. By this 
reason we work with a vector magnetic potential such that: 

B = V X A (3. 19) 

from equations (|3. 19|l and (|2. 10|l . using the vector identity 

V X (V X A) = V(V • A) - V^A (3. 20) 

we obtain the diffusion equation for the vector potential: 

47r r ^1 -> 

V^A = vx(VxA) +V(V-A). (3.21) 

L J 
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. Using the identity relation 

V{v-A) = {v-V)A+ {A- V)v + vx (V x 1) + A x (V x v). (3. 22) 

and keeping in mind that the velocity is kept constant, expression 22|l is transformed into 

vx{VxA)^V{v-A)-{v-V)A. (3.23) 
. By substituting equation H3. 23() into equation (|3. 21() we obtain 

dir r ^ ^1 -< 

V^A^ — -fia V(t7-A)-(t7-V)v4 + V(V • A). (3.24) 
L J 

The constants that appear in the first term on the right hand side of equation (|3. 24|) correspond to 
the inverse of magnetic viscosity ('7^^). Then, the equation to solve is 

r]^\/^A^\/{r],n\7 ■ A-v- A) + {v-V)A. (3. 25) 

As the velocity is constant along the x axis this equation can be written as 

7]mV^A- (t7- V)l== V(77™V • A- voA^). (3. 26) 

To solve equation (|3. 261) we use the condition V(77mV • A — vqAx) for the potential vector. From 
the solutions obtained we choose the one fulfilling the above condition for the potential vector. So 
the equation to resolve for each component of the vector potential is: 

rj^y^A^- {v-V)A^ =0 

r]^W^Ay~{v-W)Ay^O (3.27) 
77™V'A, -(z/-V)A, =0. 

It is well known that the velocity is proportional to the negative gradient of the potential function 
(py. Then we can write the last three equations as; 

7?™V^A,; - ■ \/Ai = 0, con x,y,z (3. 28) 

hence the general equation to solve is: 

?7™V^B - V(/)„ • = (3. 29) 

where we call B each one of the components x, y and, z of the vector potential. To solve this problem 
we propose a solution of the form: 

B = 0™Bo. (3. 30) 

Placing this solution into equation (|3. 29|l for the function B, we have that the potential functions 
must satisfy: 

V(/>™ = -|^V(/.,. (3.31) 

and 

V^Bo - k'^Bo = 0. (3. 32) 
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where 

fc^_!^L = 27r^^ (3.33) 

The solution of equation (|3. 321 is known as the modified Helmholtz's equation, which is given in 
terms of the spherical harmonics and the modified spherical Bessel polynomials, this is: 

oc I 

-Bo - ^ ^ [Ci„Mkr) + Di„^ni{kr)] Yi„,{0, cf>). (3. 34) 

(=0 m=-; 

As we want a converging solution, we have to require that Cim = 0. Then, the solution for Bq is: 

oo / 

Bo^Y^Y. [Di,nni{kr)]Yi,n{e,(l>). (3. 35) 

The solution of equation H3. 3f|l can be obtained by direct integration, equalizing component by 
component in both sides of the equation, we obtain: 

^„ = Cie-^'*'- = Cie^\ (3. 36) 

The integration constant Ci is included into constant Dim- Accordingly to equation H3. 3U|) we build 
a solution for the components of the potential vector A — (pmB: 

oo I 



A, = e*'-^ J2 Y [Dimniikr]] Yi^iO, 0). (3. 37) 

oo / 

^y^e'^^E E [EirrMkr)]Y,m{e,cj,). (3. 38) 



l=Q m=-l 



and 

oc I 

A, = e'^^ ^ Y [Fimuiikr)] Fz„(0, cf,). (3. 39) 

1=0 m=-l 

As in the proximities of the source the potential vector obtained must approach to the one produced 
by a dipole, it can be calculated applying Laplace's equation: 

^'^'^ ^ ^72" ^i'^(^) sm(0) = — -3- 

Ayd = ^ sin(6i) cos(?!)) = ^ (3. 40) 

A,d = 0. 

Where m' is the dipole moment. Studying the angular behavior of the solutions H3. 37|l . H3. 38|l and 
l|3. 39|l we can consider a solution of the form: 

.1 -k(r-x) myil + kr) 

~ ^ k^/^r^ 
A' = 0. 
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Where m is not a dipole moment, as it does not have the right units, but it has a dependence with 
the Earth's dipole moment Of the form m cx "^"^'gf ^^3/2°'° ■ 

However, in order to have a vahd solution, we should require the divergence of A to be given by: 



T^mV • A - VqAx = 0. 



(3. 42) 



This is the condition that we have been working. Therefore, it is necessary to add an extra term, 
coming from introducing the sum of equation 13. 38(1 into Ay, and thus the final solution for the 
magnetic potential vector is: 



Ax 
A, 



^ 3 / 2 

x{l + kr) 1 



Vkr 



-k{r — x) 



(3. 43) 



A, = 0. 

From equation H3. we can determine the magnetic field which, for our case, is: 
Q-Hr-x)^^ (3xr + k (3xr^ + r^) + P {xr^ — r^)) 



By - 



^3/2^6 

e-Hr-x)^y^ (3fcr2 + r (3 + k'^r'^)) 

g-fe(r-x)^ ((3^2 -r^)r + kr^ {3z^ - r^)) e-'^'^'-^^mk^z^r^ 



(3. 44) 



It can be shown by means of direct calculation that this magnetic field satisfies Maxwell's equations 
and the diffusion equation ((2. f2|l . Nevertheless, from the analytic expression is difficult to see the 
form that the field possesses, therefore it is convenient to observe a graphic of the magnetic field in 
the space. For doing this we should keep in mind that constant k contains the information of the 
solar wind, i.e. velocity and conductivity. 

To see the influence of velocity of solar wind upon the magnetic fleld two graphics in the plane 
y = 0, m = Id are shown: figure la shows the case of A: = 0, 005 and figure lb shows fc = 0, 01 
where u is an adequate unit of length and d is an adequate unit for m. 



u u u u I n 



1 11 11 11 1 11 1 



UJ.l.tJ.I H U 



Figure No. 1. Terrestrial magnetic field modified by the presence 0/ the solar wind; the solar wind 
comes from left to right. We have taken (a) k = 0, 005 u^^ and m — Id. (h) fc = 0, 02 and 

m — Id. 

Figures No. la and No. lb show that the terrestrial magnetic field is dragged by the solar wind in 
the direction of its movement, compressing the field lines in the face of earth directed towards the 
Sun. 
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In the other hand, as the magnetic field in the magnetosphere in known, the current of charged 
particles of the solar wind traveling by the magnetosphere can be calculated using Ampere's law 

J= —V X B. (3. 45) 

47r 

and the current is 



Jy 



27r^/kr 



CTOwze'^f^-'') (3 + 3kr + k'^r'^) 

Figures No. 2a and No. 2b show the current in 2D, using vectors to visualize its direction, just 
as we did with the magnetic field. It is clear from (|2. 81) that the current is perpendicular to the 
velocity of the solar wind, this means that it does not possess a component in the i direction, by 
this reason we draw the components y and z. The values of k and m are the same as those for figure 
No. la; m = Id, fc = 0, 005 u'^. 



Figure No. 2. Induced current into the solar wind by the presence of the terrestrial magnetic field. 
The solar wind comes out of the page. We have taken (a) In the plane x — 0, k ~ 0, 005 u^^ and 
m = Id. (h) In the plane x = 500, fc = 0, 005 u^^ and m = Id. 



In equation \2. 13() we introduce the differential identity: 

(3. 47) 

Given that all temporary partial derivatives fade for short times, equation (|2. 13|l can be written as 



dt=^''-^^^dt 



p(v ■ V)t/ = -VP + - j X S + iqS/^v. 
c 



(3. 48) 



The equation of movement (|3. 48|l . according to our Process, should not be satisfied due to the 
fact that it was Taken a constant velocity in x direction for the solar wind. We must remind that 
equation (|3. 48|l tells us about the velocity of the fluid in all space, while that of our model is valid 
just far away of the planet. By this reason, in the equation of movement (|3. 48|) we cannot consider 
a constant velocity, since it would be a meaningless equation. It is convenient to consider small 
variations in the velocities, small enough to preserve the condition that far away from the planet 
velocities are almost constant. Later we will see that this consideration is correct as the magnetic 
force vanishes quickly if we move far away from Earth. It is difficult to find a solution of equation 
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48|l . however we can center in the form of the force that the terrestrial magnetic field has upon 
the solar wind, for this we will consider the Lorentz's force: 

F ^-jx B. (3. 49) 

c 

Substituting in equation H3. 491) the values of the magnetic field and of the current induced upon 
the solar wind, we obtain the force of the magnetic field of Earth upon the solar wind. In figure 
No. 3 we show a vector diagram of the forces in the plane y = 0. In this figure the effect of the 
magnetic field upon the fiuid can be observed, trying to move it toward the "poles". In the x axis, 
which is the direction of the incident solar wind, the field tends to turn back the fiuid, this force 
generates a strong pressure that, in a better model, would modify the velocity of the fiuid making 
it to go around the magnetosphere, protecting the surface of the planet from high velocity ionized 
gases. In the same direction, the x axis, when the fiuid has passed the planet it tends to brake but 
the braking force is not strong enough. 




Figure No. 3. Force of the terrestrial magnetic field upon the solar wind in the plane y — 0. The 
solar wind comes from left to right. We have taken k = 0,005 u^^ and m — 1. 

To see the behavior of this forces as a function of the distance to the source, figure No. 4 shows 
the magnitude of forces, using the same value of distance as in figure No. la. From this figure 
results that the force between the source and the fiuid to those distances vanishes. Figures No. 
la and No. lb have an interval in x and in y of (— 500u, 500u), while the magnitude of the force 
is obtained to be relevant in the distance interval (—50m, 50u), still inside the planet as the circle 
representing Earth possesses a radius of 80u. The force has been taken in adequate units of /. 

Given the great velocity of incoming solar wind, far away from Earth it does not experience any 
deviation. 
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Figure No. 4- Magnitude of the force between the terrestrial magnetic field and the solar wind in 
the plane x = 20m. The solar wind goes downwards. We have taken k — 0, 005 u^^ and m ~ Id. 



4 Conclusions 

We have shown that the terrestrial magnetosphere is affected by the solar wind, suffering a length- 
ening of the magnetic field lines. On the other hand, due to the presence of the terrestrial magnetic 
field, certain currents arise inside the solar wind, generating an additional magnetic field which 
superposes to the dipole magnetic field, creating the total field in the magnetosphere. 

We have shown in this model that the force exerted by the terrestrial magnetic field upon the 
solar wind vanishes very quickly with distance, when we move away from Earth. As a consequence, 
our model is valid for determinate the magnetic field far from the Earth. 
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